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Strong influence of the complex bandstructure on the tunneling electroresistance: 
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The tunneling electroresistance (TER) for ferroelectric tunnel junctions (FTJs) with BaTi03 
(BTO) and PbTi03 (PTO) barriers is calculated by combining the microscopic electronic structure 
of the barrier material with a macroscopic model for the electrostatic potential which is caused 
by the ferroelectric polarization. The TER ratio is investigated in dependence on the intrinsic 
polarization, the chemical potential, and the screening properties of the electrodes. A change of sign 
in the TER ratio is obtained for both barrier materials in dependence on the chemical potential. 
The inverse imaginary Fermi velocity describes the microscopic origin of this effect; it qualitatively 
reflects the variation and the sign reversal of the TER. The quantity of the imaginary Fermi velocity 
allows to obtain detailed information on the transport properties of FTJs by analyzing the complex 
bandstructure of the barrier material. 

PACS numbers: 31.15.A, 73.23. Ad, 85.50.Gk, 73.40.Gk,77.80.Fm, 71.20.Ps 
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While in the last two decades the field of spintronics 
dominated the area of high-technology memory and sen- 
sor devices, non- volatile ferroelectric memories gain more 
and more attention 1 . Recently, two phenomena which are 
known in science for more than 80 years, namely ferro- 
electricity and quantum mechanical tunneling, are com- 
bined in ferroelectric tunnel junctions (FTJs) and could 
soon compete with common spintronic devices 2 . 

In one of the earliest concepts introduced by Esaki et 
al. , that is, a 'polar switch', the tunneling barrier was 
composed of a ferroelectric. Its spontaneous, nonvolatile 
polarization can be switched by an applied electric field. 
At the best the switching of the barrier's polarization 
would create two different electrostatic potentials and 
,therefore, could lead to two distinct levels of the tunnel- 
ing conductance. This tunneling electroresistance (TER) 
was predicted by Zhuravlcv et a/. 4,5 and could be em- 
ployed in binary-logic devices. 

While for a long time ferroelectricity was limited to di- 
mensions above several hundred nanometers and viewed 
as a collective phenomenon 6 , its combination with quan- 
tum mechanical tunneling was impossible. Hence, a 
TER effect could not be established. However, during 
the past few years, there has been tremendous progress 
in understanding finite-size effects in ultrathin ferroelec- 
tric films. Today, theory predicts 7-9 and experiment 
demonstrates 10,11 the presence of ferroelectricity for film 
thicknesses down to a few unit cells; so, both support 
the concept of FTJs. After the breakthrough of Contr- 
eras et al. 12 , a number of recent experiments 13 success- 
fully investigated the transport properties of FTJs using 
piezoresponse and conductive atomic force microscopy 
and found evidence for a TER effect. 

The purpose of this paper is to investigate theoreti- 
cally fundamental relations between the electron trans- 
port through mctal/ferroelectric/metal FTJ devices and 



the polarization-state-dependent electronic structure of 
the barrier. 

Here, we consider both ultrathin ferroelectric BaTi03 
(BTO) or PbTi0 3 (PTO) as barriers. Recently, cal- 
culations based on free-electron models 4,5,18 ' 19 or first- 
principles methods 20 explained the basics of the TER 
effect. 

In contrast to those investigations, our calculations 
successfully combine an analytical and phenomenologi- 
cal model with ab-initio calculations. This ansatz en- 
ables a systematic study of the TER in dependence on 
the barriers electronic structure and the position of the 
systems chemical potential Ep separated from interface 
effects and the influence of the metallic electrodes. 

For the electronic structure we applied multiple- 
scattering theory, namely Korringa-Kohn-Rostoker 
(KKR) methods 21,22 , to obtain an accurate state-of-the- 
art description of the electronic states (cf. Section II B). 
In all steps, we applied the local density approximation 
(LDA) to density- functional theory (DFT). The calcu- 
lation of the transport properties was realized within 
a Landaucr-Biittikcr picture 23 , involving the ab-initio 
complex bandstructure and analytical expressions for 
the electrostatic potential which is caused by the 
ferroelectric polarization and the different screening of 
the interface charges in the metallic electrodes 4,5 . This 
approach provides a flexible and fast computational 
means to observe the influence of material-specific 
parameters, especially the electrodes' properties, on 
the TER effect, while even taking into account the 
exact electronic structure of the barrier. Therefor, the 
proposed approach describes the transport quantities 
of FTJ's in a reliable way, which does not require a 
full self-consistent calculation of the whole junctions 
electronic structure. 

The introduced term of the inverse imaginary Fermi 
velocity allows to estimate the TER effect for a given fer- 
roelectric barrier material, knowing the barrier thickness, 



polarization strength and the the ratio of the electronic 
screening lengths in the electrodes. 

This paper is organized as follows. In section II A the 
construction of the electrostatic potential in the concept 
of FTJ is introduced. Section II B discusses issues of the 
self consistent electronic structure calculations, while sec- 
tion II C summarizes the calculation of the conductance. 
Section III presents the obtained results, starting with a 
detailed discussion of the complex bandstructure of BTO 
and PTO in III A and focusing afterwards on the be- 
havior of the TER in dependence of the position of the 
chemical potential (Sec. IIIB1) and the electrodes elec- 
tronic screening lengths (Sec. IIIB3). The term of the 
inverse imaginary Fermi velocity is introduced in section 
III B 2 to provide a comprehensive relation between the 
electronic structure and the obtained transport proper- 
ties. 



II. METHODOLOGY 

The electrostatic potential in the barrier was derived 
from interface charges caused by the ferroelectric polar- 
ization of the material and the electronic screening in 
the electrodes. The electronic structure and so the de- 
cay rates of the states in the electronic band gap were 
calculated sclf-consistcntly. A perfect 2D periodicity was 
assumed perpendicular to the transport direction. The 
transmission probability was calculated using a WKB 
approximation. The matching of the wavefunctions at 
the interfaces was assumed to be equal for all states, so 
it does not enter the TER ratio in this approximation. 
Interface resonances which might rarely appear in real 
junctions are not considered and do not change the re- 
sults qualitatively. 



A. Ferroelectric tunnel barrier 

The physical mechanism behind the TER effect is the 
change of the electrostatic potential caused by the rever- 
sal of the electric polarization P in the ferroelectric 4 ' 5 . 
For sufficiently thin ferroelectric films (a few unit cells 
thick), the charges at the ferroelectric/electrode inter- 
faces are not completely screened by the metallic elec- 
trodes. As a consequence, the depolarizing electric field 
E is nonzero 7,9,24 , and its sign depends on the direction 
of P. 

Due to different screening properties of the electrodes, 
e. g. asymmetric interface termination 20 ' 25 or different 
electrode materials 12,15 , the electric polarization leads to 
an asymmetry in the potential profile. Thus, the effec- 
tive potential seen by electrons changes with reversal of 
P. The averaged height of the barrier potential differs by 
AV, which leads eventually to the TER effect (Fig. 1). 

Assuming a Thomas-Fermi model of screening and a 
constant electric field in the ferroelectric barrier, the 
static electric potential profile in the electrodes and 
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FIG. 1: (color online) Schematics of the tunnel barrier poten- 
tial along the current direction for the ferroelectric insulator, 
with the barriers polarization pointing to the right (a, P+) or 
to left (b, P-) interface. The potential is a superposition of 
the electrostatic potential (colored line), the electronic poten- 
tial which determines the bottom of the bands in the two elec- 
trodes, and the potential barrier created by the ferroelectric 
insulator. A and A are the screening lengths of the electrodes, 
cf. Eq (1). 
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with z the coordinate along the current direction 4,5,7 . 
With the screening lengths A and A of the left and the 
right electrode, respectively, the screening charge per 
unit area is given by as = dp /c(A+\)+d, where P is the ab- 
solute polarization and e the relative static permittivity 
of the barrier material. The latter was chosen for BTO 
and PTO 1500 and 250, respectively 26 . For simplicity, A 
is assumed to be larger than A. The average height of 
the potential barrier depends on the orientation of the 
polarization, this is, pointing to the left (P_) or to the 
right (P+) interface. This is due to the fact that the ad- 
ditional electrostatic potential has an averaged value of 
either +\v(°)+v( d )\/2 or — \<pW+<p(d)\/ 2 . 



B. Electronic-structure calculations 

For our computational approach a sequential multi- 
code treatment was applied. Various quantities were 
carefully cross-checked among the computer codes to ob- 
tain consistent results. Reliability is achieved by numer- 
ous convergence tests. 

Firstly, we determined the fully relaxed atomic po- 
sitions according to Fechner et al. 27 within a pseu- 
dopotential and plane- wave-basis scheme using VASP 28 . 
With the ferroelectric polarization pointing along the -re- 
direction, the values of P for polar bulk BTO and PTO 
were calculated by a Berry-phase approach 29 . 

In a second step, the electronic structures of the 
considered systems, with atomic positions from the 
first step, were calculated self-consistcntly within the 



framework of density- functional theory (DFT), using 
a scalar-relativistic screened Korringa-Kohn-Rostoker 
(KKR) Green function method 21 . The spherical poten- 
tials were treated in the atomic sphere approximation 
(ASA), using the local density approximation (LDA) for 
the exchange-correlation potential 30,31 . 

In a third step, the complex band structure k z (E, fen) 
was computed within a layer KKR, using the same poten- 
tials as in the second stage 22,32 . The imaginary part of 
k z (E,k\\) will be denoted as decay parameter k(E, fey). 
As will become clear in what follows (III A), evanescent 
states play an important role in the tunneling regime. 

Applying the scheme explained in the upcoming sub- 
section, we are now able to compute the transmission 
probability which determines the conductance and there- 
for the TER properties. 



C. Conductance calculations 

In order to predict the conductance change associ- 
ated with a polarization switching, we assume so thin 
a barrier that the dominant transport mechanism is 
quantum-mechanical tunneling and the ferroelectricity is 
preserved 8,10 . Therefore, four unitcells (d w 1.6 nm) of 
BTO or PTO are considered as electric switchable bar- 
rier. 

Assuming clastic and coherent transport, the zero-bias 
conductance per unit cell area is given by the Landauer- 
Biittikcr formula 23 . 



G= 2 4T, T (E F ,k l{ ), 



(2) 



where 



a 

T(E F ,k\\) = T exp[-2 fdzK(E F -<p(z),k\\)] 



(3) 



is the transmission probability in the WKB approxima- 
tion at the chemical potential Ep . k is the smallest imag- 
inary part of the complex wavevector (cf. Section III A) 
and is calculated by means of the ab-initio complex band 
structure. The WKB approximation takes into account 
the shape of the barrier potential ip(z) and the decay 
properties of the electronic states in the barrier gap. The 
fen integration is over the whole two-dimensional Bril- 
louin zone, using at least 6000 special points in the ir- 
reducible part of the Brillouin zone. The transmission 
prefactor To comprises the influence of the wavefunction 
matching at the interfaces to the tunneling probability. 
We used a quite rough approximation, taken Tq to be con- 
stant for all energies and inplane wavevectors fen. This 
factor drops out in the determination of the TER ac- 
cording to Eq. (4). The change of the potential in the 
electrodes in the vicinity of the interfaces is not included, 



because metals with a constant density of states are as- 
sumed. 

A main feature of transport through FTJs is a conduc- 
tance asymmetry, defined as 



TER = 



G(P-)-G(P+) 

G(P-) + G(p+y 



(4) 



where G(P + ) (G(P-)) is the conductance for the polar- 
ization pointing to the right (left) interface. 



III. RESULTS AND DISCUSSION 
A. Complex band structures 

To prepare material-specific input for the transport 
calculations, it is useful to recall the importance of the 
complex band structure for tunneling through an (ferro- 
electric) insulator (see 33 for magnetic tunnel junctions). 
The imaginary part k of the complex wavevector domi- 
nates the transport properties; its product with the bar- 
rier thickness defines the exponential decay of the wave- 
functions within the barrier and determines essentially 
the transmission probability. 

In a periodic crystal, the wavevectors fc of Bloch states 
are necessarily real. Heine 34 demonstrated that a surface 
or interface state can be obtained by matching eigen- 
states of both sides of the boundary plane. This pro- 
cedure requires in addition states with an imaginary 
or — more general — a complex wave vector, i. c. so-called 
evanescent Bloch states. The wave-function matching 
across the boundary requires eigenstates with complex 
wave numbers k z at a given real energy have to be taken 
into account. 

At an interface between a metal and an insula- 
tor, prominent states are the metal-induced gap states 
(MIGS), which are itinerant in the metal electrodes but 
decay exponentially into the insulating barrier. The dis- 
persion relation of the evanescent states is the so-called 
complex band structure. For a planar tunnel junction, 
the periodicity along the interfaces (barrier) requires that 
the projection k» of the wave vector k onto the interface 
is real and conserved. However, the perpendicular com- 
ponent fc, +in can be complex and depends on the region 
(electrodes, barrier). 

From Figure 2 we deduce fundamental band gaps of 
ferroelectric BTO and PTO of 1.79 eV and 2.89 eV, re- 
spectively. Both are increased, by 4.1% for BTO and 
21.4% for PTO, as compared to the calculated paraelec- 
tric case. The strong localization of the A 2 ' states pre- 
vents larger band gaps. 

The para-to-ferroelectric phase transition lowers the 
symmetry from cubic (Pm3m) to tetragonal (P4mm). At 
r, the bands can be decomposed with respect to the irre- 
ducible representation of the point group 4mm, as shown 
in Fig. 2. The t 2g bands split to form a doubly degener- 
ate A5 band (d ZXl d zy ) and a A2 band (d xy ). Similarly, 
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FIG. 2: (color online) Symmetry-resolved complex band 
structures for ferroelectric (a) BaTiC>3 and (b) PbTi03 at 
r in the 2D Brillouin zone (fey = 0). Middle panels show the 
conventional bands structure (Im k z = 0) along the T — Z- 
line. Left and right panels display imaginary bands of the 
first (k z = 0) and the second kind (k z — — ), respectively. 



the e g bands split into a Ai band (d z 2) and a A2' band 

(d x 2_ y 2). 

According to Chang 35 , complex bands can be classified 
as follows, (i) Real bands correspond to the conventional 
band structure and have Im k z = 0. Thus, the wave func- 
tions are the Bloch states, (ii) Imaginary bands of the 
first kind have Re k z = and Im k z ^ 0. (Hi) Imaginary 
bands of the second kind have Re k z = */c and Im fc 2 ^ 0. 
(iv) Complex bands have Re k z ^ , Re k z ^ ff /c and 
Im k z ^ with c being the periodicity in z-dircction. 

At r, the three states with the smallest decay rates 
within the fundamental band gap consist of a A 5 doublet 
and a Ai singlet for both perovskite ferroelectrics. Nev- 
crthlcss, we address, that the smallest decay rate for BTO 
is formed by either the A5 or the Ai imaginary band 
of first kind (cf. 2(a), band crossing in the left panel), 
while for PTO the smallest decay parameter is always 
associated with a Ai-like imaginary band of the second 
kind (cf. 2(b), right panel). Importantly, these complex 
bands — and hence the electric conductance — are highly 
sensitive to the ferroelectric displacements; the latter, in 
turn, depend significantly on the polarization direction 
and on the potential profile in the barrier. 

Since at each fcii the transmittancc is determined by the 
least decaying wavefunction, we address the fc||-rcsolved 
smallest imaginary part k of the wave vector at different 
energies in Fig. 3. For BTO, the main contribution to 
the transmission is expected to come from a small an- 
nulus around T for all energies within the gap. This is 
caused by the smallest decay rates near the Brillouin zone 
center. Our findings corroborate those of Velev et al. 36 . 




FIG. 3: (color online) Smallest imaginary part of the complex 
wavevector of BaTiOa (upper row, a-c) and PbTiOs (lower 
row, d-f) in the 2D Brillouin zone (edges of the Brillouin 
zone at ±""/ a )- The imaginary parts are shown for different 
energies relative to the valence band maximum: (a) 0.77 eV, 
(b) 1.30 eV, (c) 1.77 eV, (d) 1.10 eV, (e) 1.40 eV, (f) 2.10 eV. 



In contrast to the paraelectric phase (not shown), the 
smallest imaginary part of the wavevector is not located 
exactly at the Brillouin zone center but slightly off-set 
for energies close to the valence band edge (Fig. 3(a)). 
This is explained by the lifted degeneracy of the Ai and 
A5 bands at the T-point. The latter supports the im- 
portance of a wave- vector analysis in the entire Brillouin 
zone. 

For PTO, the situation is more complex in comparison 
to BTO. In particular, the smallest k are located at F 
and X edge centers of the 2BZ for energies close to the 
valence band maximum and the conduction band mini- 
mum (3(d)-(f)). In the middle of the gap, states from the 
Brillouin zone center dominate the transport (3(e)). In 
contrast to BTO, one depicts for PTO nearly the same 
damping around M (corners of the Brillouin zone), while 
sizably larger k in the rest of the Brillouin zone can be 
stated (note the different scale in top and bottom pan- 
els of Fig. 3). Thus, one expects remarkably different 
fen -transmissions and smaller tunneling currents as com- 
pared to BTO. 



B. Transport Properties 

Having provided complex band structures, the tunnel- 
ing conductance of the FTJs is investigated in depen- 
dence on the polarization of the barrier, the screening 
lengths in the electrodes, and the chemical potential with 
respect to the fundamental band gap of the barrier. For 
ferroelectrics, the Fermi energy can be adjusted within 
the gap by crystal doping 37 , while a shift of the chemical 
potential can be realized by a small bias. 



1. Energy dependence of the TER 

We recall that the slope of the imaginary part of the 
complex wavevector k and the K distribution in the Bril- 
louin zone change remarkably with energy (cf. Figs. 2 and 
3). As a consequence, we focus on three positions in the 
gap: close to the valence band, in the center of the gap 
and close to the conduction band. 

In Figure 4 the calculated conductance and TER for 
FTJ's with BTO (upper row, a-c) and PTO (lower row, 
d-f ) as a ferroelectric barrier material are presented. In 
the upper part of each panel the conductance referring to 
Eq. (2) is shown. Full blue (open red) circles denote the 
dependence of the conductance for the polarization point- 
ing to the left (right) interface, while the ratio of screen- 
ing lengths of the electrodes are fixed at A /x = 9. This 
might refer to a noble metal electrode (AA g = 0.07 nm) 
and a ferromagnetic electrode (AsrRu0 3 = 0.6 nm) or a 
half metallic electrode (AL ai _ x Sr x Mn0 3 = 0.2 — 1.9 nm) as 
reported in 4,7 ' 13 . In the lower parts of the panels 4 (a)- 
(f) the TER determined by Eq. (4) in dependence on the 
abolute value of the polarization is represented by open 
black circles. 

As we consider different energies in the gap — close to 
the valence band (a,d), near the gap center (b,e), and 
close to the conduction band (c,f) — one easily recognizes 
the increase of the conductance near the gap edges. This 
is readily explained by the imaginary part of the wavevec- 
tors which decreases for certain fey (compare Figs. 2 and 
3). Hence, the transmission probability and the con- 
ductance (ref. to Eq. (3)) increase. In a free electron 
picture as used in previous explanations of the TER 
effect 4,5,18,19 , the conductance would decrease continu- 
ously by shifting the chemical potential from the band 
edges. Our approach takes into account a more realistic 
electronic structure of the ferroelectric barrier. 

The behavior of the conductance for different chemical 
potentials can be understood qualitatively by the com- 
plex band structure at T (see Fig. 2). For energies near 
the valence band the P + -conductance is larger than the 
P conductance (compare Fig. 4(a,c)). This is caused 
by the lower decay rates for the electron wave functions 
contributing to the tunneling current in the P+ state. 
The shift of the bands to higher energies due to the pos- 
itive contribution of the electrostatic potential leads to 
smaller values of the imaginary part of the wave vector 
at the chemical potential. Taking into account the defini- 
tion of the TER (ref. Eq. (4) a positive TER occurs. For 
energies close to the conduction band edge the behavior 
of the conductances and, thus, of the TER reverse, as 
does the curvature of the imaginary bands reverse, too. 

For bulk polarizations of Pbto = 22.9/iC/cm 2 and 
Pp TO = 94.3/zC/cm 227 we achieve high TER values 
of nearly 3% (E F = E v + 0.77 eV) and 27% {Ep = 
Ey + 1.10 eV) for chemical potentials near to the va- 
lence band and allmost 6% (E F = E v + 1.77 eV) and 
61% (Ep = E v + 2.10 eV) for energies close to the 
conduction band. The value for PTO is in agreement 



with experimental values of nearly 60 % for a FTJ with 
Pb(Zr,Ti)0 as ferroelectric barrier with different elec- 
trode materials 12 . Recent findings for FTJs based on 
strained BTO show a TER of 16 % 15 , which is in aston- 
ishing agreement with our predictions (see Fig. 4(c)) for 
strained thin films (note that then Pbto = 50/iC/cm 2 
according to 38 ). 

As for energies lying close to the gap edges the discus- 
sion can be made in a generally more intelligible pic- 
ture, the situation is quite more complex for energies 
in the middle of the gap. In figure 4(b) the evaluation 
of the TER for a FTJ consisting of BTO is shown for 
Ep = Ey + 1.30 eV near the point where the symmetry 
character of the dominating evanescent wave functions at 
the T-point changes from Ai to A5 (compare Fig. 2(a), 
left panel) and consequently the slope of the imaginary 
bands change sign. 

For polarization strengths less than Pbto = 
50 //C/cm 2 the TER nearly vanishes, while for higher po- 
larization strengths the conductance for the P state in- 
creases and a negative TER of —15% can be reached. 
As noted earlier, those higher polarization strengths 
might be achieved in stressed thin films where values of 
Pbto = 60/zC/cm 2 and even Ppto = HO/iC/cm 2 were 
reported 38 . In the case of PTO the TER is barely visible 
for the whole range of polarizations with values below 
1% for an energy of Ep = Ey + 1.40 eV. This is near 
the center of the complex wave vectors loop as shown in 
Fig. 2(b), right panel. Close to the maximum of the band 
with respect to the decay rates, a small shift of the po- 
tential caused by the electric polarization reversal causes 
a negligible change of the decay rates, only. 

As the decay parameters of all states in the Brillouin 
zone are comprised for the calculation of the transport 
properties, different curvatures of k are involved in the 
total resulting TER and a straightforward explanation of 
the properties of the TER is hardly possible. In the pre- 
ceding, we established major trends by discussing three 
exemplary energies. Now, we turn to the TER in depen- 
dence on the position of the chemical potential relative 
to the barrier gap. In the upper row of Fig. 5, the most 
remarkable feature is a change of sign of the TER at 
energies of about 1.5 eV for both materials. This result 
cannot be explained by a free electron model, where the 
negative slope of the imaginary band would lead to neg- 
ative values of the TER over the whole range of the band 
gap. Furthermore, nearly free electron models or tight- 
binding approaches would produce only one sign change 
at an energy approximately in the middle of the gap, 
which is at about 0.8 eV and 1.4 eV for BTO and PTO, 
respectively. 

For the chosen polarization values and for energies 
close to the valence band of BTO, the TER is expected 
to be as twice as large as for energies close to the con- 
duction band. For PTO even higher values of 40 % (close 
to the valence band) and -60 % (close to the conduction 
band) arc expected. 



o 



a 

7 
6 
S 
4 
3 

10 



t 

H -20 
-30 

-? a 

= 

J* 7.5 

'o 

C 6 

o 

5 3 

1 

LJ 1.5 

PC 

| ° 
a -20 

" -40 

RC 




(a) 



I ■ ■ I ' ■ I ■ ■ I ' ■ I 







*(d) 




i- 




I I I I " j I I I I I I I I [ 



I ■ ■ I ' I I ' ' I ' ' i 




aoooooooooopo 



» 




■I' M 



I' M 





' ' I ' ' I ■ ' I ' ■ I ■ ' I . 




20 40 60 50 100 20 40 60 80 100 20 40 60 60 100 



Polarization (nCtoir) 



Polarization (^C/cm") 



Polarization (jiOcrrO 



FIG. 4: (color online) Conductances and tunneling electroresistances (TERs) for the two ferroelectric configurations of a FTJ 
with BaTiOs (upper row, a-c) and PbTiC>3 (lower row, d-f) as a ferroelectric barrier. The barrier thickness was fixed to 
four unit cells (~ 1.6 nm). Both, conductance and TER, are shown in dependence on the absolute value of the electrical 
polarization as well as for different chemical potentials with respect to the valence band maximum: (a) 0.77 eV, (b) 1.30 eV, 
(c) 1.77 eV, (d) 1.10 eV, (e) 1.40 eV, (f) 2.10 eV (that is, as in Fig. 3). Full blue (open red) circles denote conductances for 
the polarization pointing to the left (right) interface. The inset in (e) depicts the conductances in a smaller scale, showing a 
small but pronounced variation. A and A are the screening lengths in the left and right electrode, respectively , with fixed ratio 
A A = 9. 



2. Inverse imaginary Fermi Velocity 

A qualitative explanation for the previous described 
energy dependent behavior of the TER is provided by 
the introduction of the inverse imaginary Fermi velocity 



v F \E) = h 



, dn.(k h E) 
1 dE 



BZ- 



(5) 



This quantity can be seen as the fcn -averaged energy 
derivative of the decay parameter k of the complex band 
structure in the Brillouin zone. If this parameter is small, 
the additional electrostatic potential caused by the elec- 
trical polarization P_ or P + does not change the effective 
decay parameter remarkably. So, a small TER effect is 
expected. On the contrary, a large energy derivative of 
the decay parameter k which most likely occurs close to 
the band edges points to a large TER effect. 

To illustrate the correlation of the inverse imaginary 
Fermi velocity and the TER, an approximation for a sin- 



gle feii is given in the following. Small polarizations will 
lead to small differences in the complex bandstructure for 
both P±. Hence, the imaginary part «_ of the complex 
wavevector k z for P_ can be written as 



dE 



(6) 



by expanding the imaginary part «+ for P+. The TER 
given in Eq. (4), including Eqs. (3) and (6) is directly 
heading to the approximation 



Here. 



TER«a- (h-vp)- 1 ■ P-d 2 



e(A-X) 



e [e(A + X) + d]' 



(7) 



(8) 



is a unique constant for each junction, as e and d contain 
material-specific information of the ferroelectric barrier 
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FIG. 5: (color online) Tunneling electroresistance (TER) ver- 
sus chemical potential for BTO (a) and PTO (b) relative to 
the valence band maximum of the ferroelectric. The screening 
lengths' ratio is fixed to 9. In the lower part the correspond- 
ing inverse imaginary Fermi velocity for BTO (c) and PTO 
(d) is depicted. For further details see text. 



and the screening lengths A and A are determined by the 
electrodes. An explicit value of a for BTO (A = 6 nm, 
A = 0.7 nm, e = 1500, d = 1.6 nm) would be a « 60 • 
10 6 eVm/C. 

The effect of all states in the whole Brillouin zone has 
to be recognized to describe the TER effect in a reliable 
way. Therefor the inverse imaginary Fermi velocity is a 
combined quantity for the damping of all wavefunctions 
in the Brillouin zone and hence a value for the probability 
of transmission. The lower panels in Fig. 5 depict the in- 
verse imaginary Fermi velocity for ferroelectric BTO (b) 
and PTO (d). The behavior of the TER is fully repro- 
duced qualitatively and quantitatively. To be more spe- 
cific, the energy position of the sign change in the TER 
as well as the magnitude of the TER for different chem- 
ical potentials can be achieved by analyzing the inverse 
imaginary Fermi velocity. The different ratios between 
the TER and the inverse imaginary Fermi velocity are 
caused by the relative static permittivities for BTO and 
PTO. This also prevents the calculation of the TER for 
PTO junctions for chemical potentials close to the band 
edges, as electrical breakthrough and metallic behavior 
would be the consequences. 

The perfect agreement of the behavior of TER and 
the inverse imaginary Fermi velocity is to some extent 
accidental. The average over Vp is taken in the whole 
2BZ, whereas the transmission defined in Eq. (3) is dom- 
inated by the fey states with the lowest decay parame- 
ters k(E, fey). It seems that these regions of the 2BZ 
are also dominated by the averaged inverse Fermi veloc- 
ity. For the barrier materials BTO and PTO the latter 
could be found. Introducing the exact material proper- 
ties of the electrodes, the symmetry selection caused by 
the wavefunction matching conditions at the interfaces 
would most likely weaken the link between the TER and 
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FIG. 6: (color online) TER dependence on the polarization 
strength and the position of the chemical potential for FTJ 
with (a) BTO or (b) PTO barriers. The screening lengths' 
ratio was fixed to a /a = 9. 



To complete the discussion, we show in Fig. 6 the de- 
pendence of the TER on P and Ep. For BTO, the en- 
ergetic position of the sign change of the TER varies 
slightly with the polarization strength, which is most 
likely linked to the change of symmetry character of the 
dominating evanescent wave functions associated with a 
sign reversal of the inverse imaginary Fermi velocity at 
the crossing point of the two imaginary bands of the first 
kind at an energy of about 1.4 eV above the valence band 
edge (compare top panel in Fig. 2). In the case of PTO 
the sign change is obtained for a chemical potential of 
about 1.4 cV above the valence band edge, independent 
on the intrinsic polarization strength. Again it is conspic- 
uous that the highest values for the TER occur near the 
band edges, where the inverse imaginary Fermi velocity 
is maximized. The largest absolute values for the TER 
are found near the valence band maximum or conduction 
band minimum for BTO or PTO, respectively. 



3. Influence of the electrodes ' screening lenghts 

To get more insight on how the TER depends on the 
electrodes' properties, we present in Fig. 7 its depen- 
dency on the ratio of the screening lengths for Ppto = 
70//C/cm 2 . The spreading of the conductances for P_ 
and P + increases with the difference of the screening 
lengths. An interchange of the electrodes will not af- 
fect the conductances, as only the polarization direction 
is reversed. An asymmetry of the screening lengths in the 
electrodes is mandatory to obtain a nonzero TER. For 
equal screening lengths, one obtains G(P-) = G(P+) 
and, hence, the TER vanishes. On the other hand, an 
interchange of the screening lengths in the electrodes re- 
verses the sign of the TER as a result of the definition 
in Eq. (4). If one interface is a free surface, the ratio of 
the screening lengths will be maximized, leading to giant 
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FIG. 7: (color online) The conductance for the two ferroelec- 
tric polarization states and the resulting TER of a FTJ with 
PbTiOa as a ferroelectric barrier material is shown. Full blue 
(open red) circles denote the dependence of the conductance 
for the polarization pointing to the left (right) interface. The 
absolute polarization is fixed to Ppto = 70 fiC/cm 2 . Both, 
conductance and TER, are shown in dependence on the ratio 
a /a of the screening lengths of the electrodes material and for 
different chemical potentials with respect to the valence band: 
(a) 1.10 eV, (b) 1.40 eV. 



values of the TER as observed with conductive atomic 
force microscopy on BTO and PTO surfaces 14,16,17 . 



dependence on the position of the chemical potential, 
the electrical polarization, and the ratio of the screening 
lengths in the electrodes. In dependence on the chem- 
ical potential position a sign reversal of the TER ratio 
is obtained. This is caused by a sign change of the en- 
ergy derivative of the decay parameter k. The latter was 
comprised by the inverse imaginary Fermi velocity. The 
behavior of this quantity in the band gap is in agreement 
with the TER ratios. 

As a consequence of the presented results we demon- 
strate that the analysis of the complex band structure is 
a powerful tool to gain useful insight into the behavior 
of the TER effect more precise than the analysis of free 
electron like models and more efficient than calculations 
of the full contact geometry done in previous works. 

As a sign change of the TER in dependence on the 
chemical potential is achieved, an experimental proof by 
application of different electrode materials or doping lev- 
els is proposed. We again emphasize that the shift of the 
chemical potential and thus the desired sign change of 
the TER might be also achieved by a small applied bias. 
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IV. CONCLUSION AND OUTLOOK 

The tunneling electroresistance in ferroelectric tunnel 
junctions with BTO and PTO barriers was discussed in 
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